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ABSTRACT 


Without  restriction  to  plane  or  inflniteslnally  small  motions,  the 
differential  equations  for  a  spinning  rotationally  symmetric  satellite 
have  been  established,  using  new  methods,  in  the  form  of  two  second- 
order  differential  equations.  An  Integral  of  these  equations  has  been 
found. 
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GYROTHEORY  OF  A  SPINNING  ROTATIONALLY 
SYMMETRIC  SATELLITE 

A  New  Integral  of  the  Equations  of  Motion 


Although  some  of  the  following  deductions  can  be  found  dispersed  In 
various  American  and  Gennan  textbooks,  the  whole  approach  given  in  this 
paper,  as  well  as  many  results,  seems  to  be  new.  It  is  worth  vdille  to 
go  back  in  the  development  to  the  very  beginnings  of  theoretical  physics, 
i.e.,  Newton's  laws  and  Hamilton's  principle. 

We  are  particularly  interested  in  the  motion  of  the  spinning  satellite 
imder  the  influence  of  the  gravitational  field  of  the  earth,  which  shall 
be  assumed  as  spherically  symmetric.  Lagrange-Hamil ton's  theory  immedi¬ 
ately  yields  the  equations  of  motion  for  N  mass  points; 


3L  d  8L 

9xj|^  dt  9x^ 


9L  .  d_  |L_ 
9y^  dt  9yj^ 

9L  _  d  9L 
dt  dz^ 


0 


0 


0 


where 


(1) 


^  “  Z  (  Z  ’'ik  (  *i  -  *k  )  (2) 

i  i,k»l 

(The  vector  3^  has  the  con^ponents  (xj^,  y^,  Zj^)) 


means  the  gravitationed  potential  of  the  mass  point.  The  constants 
are  Lagrangian  multipliers  corresponding  to  the  bovindary  conditions 

(xi  -  xjj)®  "  const. 

which  express  the  fact  that  for  a  rigid  body  the  distance  between  two 
etrbitrary  mass  points  remains  constant. 

Equations  (l)  can  now  be  condensed  into  the  vector  equation 


N 


mi 


Xi  +  Xi  +  2  y  “  0 


k=l 


We  substitute 


Xi  -  Xq  +  yi 

and  neglect  hi^er  than  linear  terms  of 


1^1 


(5) 


(6) 


which  is  certainly  Justified  because  the  diameter  of  the  satellite  is 
extremely  small  compared  with  the  radius  of  the  earth.  From  equation  (5) 


it  follows 


=  (^’o  +  *  ^o*  +  2ro  »  Xq®  +  2^  y^ 


or 


XiYi 


(7) 


Therefore,  equation  (^)  can  he  transformed  into 


N 


“i  ^  +  “i  y'i  *  -  ^  ^  yi)*o  "  ^  yi  "  2  ^  -  y^) 


k-1 


(8) 


We  now  sum  over  all  mass  points  and  take  into  account 


N 

Y  m^  =  0  (9) 

i=l 

which  states  that  3co  is  the  center  of  gravity  of  the  satellite,  and 

N 

Y_^  ^ik  (yi  -  yk^  “  0 

i,k=l 

Equation  (lO)  is  a  consequence  of  the  fact  that  since 

Xki  are  multipliers  belonging  to  the  same  boimdary  condition 

(3cj^  ~  xjj)^  =  const. 
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We  obtain 


•  • 

(E  mi)xo 


|i(E  IIL,  )  _ 

- T- 


(11) 


Equation  (11 )  states  that  the  motion  of  the  center  of  gravity  of  the 
satellite  is  not  influenced  by  the  satellite's  librations,  at  least  not 
as  long  as  we  neglect  second  and  higher  powers  of  ly^l/  1^1*  From 
equations  (8),  (9)^  and  (11 )  we  conclude 


iz  5n  m.  ^  nu  _ 

”i  yi  “  +  — ~  (^  yi)^  -  yi 

r  5  r  3 

o  o 


N 

2  ^  Mk  (yi  -  yfc) 

k»l 


(12) 


If  we  now  form  the  vector  product  of  (12)  with  y^^,  sum  over  i, 
and  take  into  account  that  we  obtain 

Z..  'T''  m.  _ 

“i  (yi  X  yj^)  =  +  2^  - ^  (Xq  yj^)  (xq  X  yj^)  (15 ) 

^o 

or 

^  (yi  X  yi)  -  +  ^  (3Eq  y^)  x  ^i)  {ik) 

We  now  make  use  of  a  fundamental  formula  describing  the  motion  of  a  rigid 
body; 

yi  “  yi  X  V  (15 ) 


k 


where  v,  the  Instantaneous  vector  of  rotation,  is  Independent  of  1. 
Substituting  (15)  in  equation  (l^^)  we  obtain 

dt  Z  [(yi^)v  “  (v  yi)yi]  ■  +  Z  ^ 

*"0 

For  a  short  time  we  shall  leave  the  vector  representation  and  shall 
introduce  the  principal  euces  of  inertia  through  the  center  of  gravity  as 
coordinate  axes.  Because  of  the  rotational  symmetry  of  the  satellite, 
two  of  these  sixes  are  determined  only  to  a  rotation  about  the  third  sucis 
corresponding  to  the  suds  of  rotation.  They  may  be  chosen  arbitrarily 
but  perpendicular  to  each  other  and  to  the  third  axis. 

Let  Vj^,  Vg,  Vg  be  the  components  of  v  with  respect  to  such  a 
coordinate  system,  1^^  being  the  moment  of  inertia  corresponding  to  the 
axis  of  rotation,  12  =  13  being  the  moment  of  inertia  of  an  sixis  througpi 
the  center  of  gravity  and  perpendicular  to  the  eucis  of  rotation.  Let  Xj^, 
x^}  be  the  con^nents  of  x^  with  respect  to  our  coordinate  system. 
Equation  (16)  can  then  be  written 

it  ° 

It  (la  ’'2)  -  -  dl  -  la)  (IT) 

o 

~r  da  V3)  -  +  Xa  (li  -  Ig) 

dt  rg® 
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We  introduce  the  vectors 


e  »  (vi,  0,  0) 

h  »  (0,  Vg,  Vg) 

and.  condense  equations  (17)  into  the  vector  equation 

^  (Ii  -  I2)  _  _  _ 

Ii  a  +  I2  h  - - ::: -  (x^  a)  X  a)  (18) 

To®  a  2 

We  still  need  addit5.onal  equations  in  order  to  eliminate  the  indeter¬ 
minacy  of  h.  At  first  we  have 

a  ii  =  0  ■  (19) 

Secondly  we  remark  that  the  end  point  of  the  unit  vector 

(20) 

whose  initial  point  coincides  with  the  center  of  gravity,  is  a  point 
rigidly  connected  with  the  gyro  satellite*  Therefore,  formula  (15)  can 
be  applied  to  e: 

e  =  eXvBeX(a+h)»eXh 


or 


a  ^ 


(21) 
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From  equations  (19)  aiid  (21)  we  conclude 


a  h  =  a  h  s  0 

and  forming  the  scalar  product  of  equation  (l8)  and  a 

a®  =  const.  » 


(22) 


(25) 


We  replace  3^  "by  rgW  where  w  is  a  unit  vector  indicating  the 
direction  earth  center  -*■  satellite. 

Combining  equations  (l8),  (20 ),  (21),  and  (25),  we  obtain  the  funda¬ 
mental  differential  equation  for  the  motion  of  the  satellite's  axis  of 
rotation,  represented  by  the  unit  vector  e: 


Ii  %e  +  I2  (e  X  e)  = 


(w  e)  (w  X  e)|  (2h) 


Equation  (24)  is  valid  for  any  elliptical  orbit. 

After  determining  e  from  equation  (24)  h  can  be  derived  from 
equation  (21)  by  differentiation  and  simple  algebraical  processes. 
Since 

V  »  cOqC  +  h 

we  master  the  gyro  motion  of  the  satellite. 

Let  us  apply  formula  (24)  to  the  case  of  a  circular  orbit  \^ere 
r^  w  const.  In  order  to  simplify  the  formulae  we  first  introduce 

T  =  fls  t 


7 


as  a  new  independent  variable  where  flg  means  the  angular  velocity  of 
the  satellite's  rotation  about  the  earth.  Formula  {2k)  assumes  the  form 


«o  «s  e'  +  ^2  ®) 

where  the  primes  denote  differentiation  with  respect  to  t.  Using  the 
abbreviations 


311  (li  -  la) _ _  _ 

i-li - ^  (we)  (wxe)  (25) 


‘"o  ^s  “ 

» 

I2  «s^  -  i  (26) 


-  1%  di  -  la)  -  K 

equation  (25)  can  be  written 

e*  +  Ig  (e”  Xe)=K(we)  (wXe)  (27) 


It  is  advantageous  to  refer  e  to  a  coordinate  system  which  rotates 
about  the  earth  in  the  same  way  as  the  satellite *s  center  of  gravity. 

The  axes  of  this  coordinate  system  may  be  defined  by  the  unit  vectors  w^ 
n  where  n  is  defined  by 


n-wx'^^^wxt 

t  «  w’ 


(28) 


The  unit  vector  n  is  a  constant  vector  representing  the  normal  vector 
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of  the  orbit  plane.  We  express  e  as  a  linear  combination  of  w,  t,  n 


e  *  >vw  +  >it  +  Pn 


since 


we  have 


w'=t,  t'»-w,  n*»0 


(29) 


e*  M  -  h)w  +  (ji*  +  'A)t  +  P*n  (30) 

If  we  move  with  the  satellite  and  observe  the  attitude  of  its  axis 
of  rotation  vdiich  is  given  by  the  unit  vector  (?v,  P)  we  observe  that 
this  vector  moves  on  the  surface  of  the  unit  sphere  at  a  velocity  the 
vector  of  \diich  has  the  components  A',  P’  in  the  satellite  system. 
Equation  (30)  shows  the  relation  between  the  velocity  vector  e'  ob- 
seived  in  a  universe-fixed  system  and  the  velocity  vector  (>v' ,  p',  P') 
observed  in  a  satellite-fixed  system  can  be  written 

e*  a  z'  +  (n  X  z)  (31) 

where 

z  =  P,  P) 

and  +  p^  +  P^  a  1. 

Applying  the  operator  of  equation  (31 )  . 

+  (S  X  I) 
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twice,  we  get 


e"  =  z"  +  2(n  X  z')  -  z  +  (z  n)n  (32) 

Equations  (31 )  and  (32)  enable  us  to  write  the  equation  of  motion 
(27)  In  the  satellite -fixed  system  w,  t,  n; 

^  z*  +  (n  X  z)  ^  +  I2  ^  (z"  X  *)  +  2(n  z)z’  +  (z  n)  (n  x  z)  ^ 

(35) 

«  K  (w  z)  (w  X  z) 

From  equation  (35  )j  apparently  more  complicated  than  equation  (27), 
an  Integral  of  the  equation  of  motion  can  easily  be  derived,  by  forming 
the  vector  product  of  (33)  and  z' .  We  obtain 

(n  z*)  +  ig  (z'  z")  +  Ig  (n  z)  (n  z')  =  K  (w  z)  (w  z’)  (3^) 


Integration  of  equation  (3^)  yields 


(n  z)  +  i  I2  ^  z'^  +  (n  z)^  y  “  2  ^  “  const.  (35) 

Let  us  now  write  the  equation  of  motion  (33)  In  the  coordinate  system 
X,  p,  P.  The  motion  can  be  described  by  three  algebraic  differential 
equations.  One  of  these  equations  is  identical  with  the  integral  (35 )^ 
the  second  equation  is  obtedned  by  forming  the  sceilar  product  of  equation 
(33)  and  n,  and  the  third  equation  expresses  the  fact  that  the  vector 
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(X.,  jx,  P)  is  a  vinit  vector; 


Ii  P  +  I  ^2  +  P'^  +  p2)  -  1  KX.®  »  const. 


P'  +  I2  (^'V  -  ^")  +  2I2  PP'  =  KX+i  (56) 


X.®  +  tx®  +  P®  -  1 


Introducing  polar  coordinates 


X.  St  cos  <|>  sin  6 


^  <3  sin  (|>  sin  6 
P  «  cos  6 


(37) 


The  system  (36)  is  equivalent  to 


ijL  cos  0  +  I  I2  (sin®  0  <j)'®  +  0'®  +  cos®  0) 


—  K  sin®  0  cos®  (b  =  const. 
2  ^ 


(38) 


(cos  0)'  -  Ig  (sin®  0  <|)' )’  +  I2  (cos®  0)'  ts  k  sin®  0  cos  «t*  sin  <|> 
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